Analytic treatment of the precessional (ballistic) contribution 
to the conventional magnetic switching 
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We consider a switching of the magnetic moment with an easy axis anisotropy from an "up" to 
a "down" direction under the influence of an external magnetic field. The driving field is applied 
parallel to the easy axis and is continuously swept from a positive to a negative value. In addition, 
a small constant perpendicular bias field is present. It is shown that while the driving field switches 
the moment in a conventional way, the perpendicular field creates an admixture of the precessional 
(ballistic) switching that speeds up the switching process. Precessional contribution produces a 
non-monotonic dependence of the switching time on the field sweep time with a minimum at a 
particular sweep time value. We derive an analytic expressions for the optimal point, and for the 
entire dependence of the switching time on the field sweep time. Our approximation is valid in a 
wide parameter range and can be used to engineer and optimize of the magnetic memory devices. 



I. INTRODUCTION 



Conventional magnetic switching by an externally ap- 
plied magnetic field H is the basis of magnetic recording 
in hard disk drives and other devices. The speed at which 
the moments of the magnetic bits can be switched be- 
tween the two easy directions has obvious implication for 
the technology performance, setting the limit for the in- 
formation writing rate. In general, the magnetic switch- 
ing time r m depends on the material parameters, sample 
size and shape. Here we are interested in the depen- 
dence of T m on the rate of change of the driving mag- 
netic field. We will call the time r/, required to flip 
the external magnetic field a "field sweep time". This 
time is of course finite in any real device used for mag- 
netic writing. Clearly, very long field sweep times will 
make the magnetic switching very slow. One can ar- 
gue then that, since T m will decrease with decreasing r^, 
the best case scenario for switching is the instantaneous 
flip of the external field with = 0. However, it was 
found numerically^ that in realistic conditions the func- 
tion T m (Th) is not monotonic and has a minimum at a 
particular value of the sweep time t£. The field sweep 
time corresponding to this minimum is optimal and any 
further decrease of will be counterproductive in terms 
of the technology performance. An analytic expression 
for r m (jh ) was announced in Ref . 0. That paper also ex- 
plained the physical reason behind the existence of the 
switching time minimum. In the present paper we pro- 
vide the detailed derivation of the approximate analytic 
expressions for the function T m (Th) and the optimal field 
sweep time r£ . We then discuss the limits of their valid- 
ity and compare analytic approximations with the exact 
numeric results. 

Before proceeding to the derivations, we would like 
to place the phenomenon under investigation into con- 
text. Magnetic switching and its speed are important 
technological parameters and were studied by many au- 
thors. We will concentrate on a single-domain magnet 



described by a moment M with an easy axis anisotropy 
energy Its switching under a static applied field is de- 
scribed by the Stoner-Wholfarth astroid^ in the f/-space 
(Fig-IU z is pointing along the easy axis). If the applied 
field is changed infinitesimally slowly, and the thermal 
fluctuations of the moment can be neglected, magnetic 
switching is achieved when the trajectory of the driving 
field in the ii-space crosses the astroid boundary The 
switching process begins at the moment of crossing and 
requires a finite time that depends on the crossing point. 

If the magnitude and the direction of the applied field 
are allowed to change during the magnetic switching pro- 
cess, the number of switching scenarios becomes infinite. 
In the most general case one wants to optimize the depen- 
dence H(t) so as to minimize the switching time. Differ- 
ent minimization schemes are discussed theoretically^— 
but their experimental realization is distant as they all 
require external fields that change exceedingly rapidly 
and have carefully controlled magnitude and direction. 

A large body of research is devoted to a restricted case 
of pulsed applied fields with fixed direction. The easi- 
est pulse shape to produce experimentally is a field jump 
from an initial value to a final value Hf. The jump 
can be instantaneous (step function) or have a certain 
rise time (smeared step function). After the jump the 
field stays at the final value H/. The equilibrium di- 
rections of the moment before and after the jump are 
described by the Stoner-Wohlfarth picture with H = H; 
and H = H/ respectively. Before the field jump the mo- 
ment resides in the minimum M; of the magnetic energy 
corresponding to the initial field. After the field jump 
the moment eventually reaches the equilibrium M/, cor- 
responding to the minimum of the energy at the final ex- 
ternal field. (We will consider situations where the energy 
has only one minimum). The switching time depends on 
both the initial and final directions of the moment or, 
equivalently, on both Hj and H/. 7 This leads to some 
interesting properties even for the idealized case of an 
instantaneous jump. For example, switching below the 
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Stoner-Wholfarth limit becomes possible^— and is ac- 
companied by interesting counter-intuitive phenomena. 
The usual property of switching time to decrease with 
with the increasing magnitude of Hf, observed for Hj 
outside of the astroid, may reverse when H/ is inside the 
astroid.— 

If instead of the held jumps one uses pulses with the 
field switched from zero to a fixed value for a finite time 
period, more options for switching time minimization be- 
come available. In the absence of magnetic anisotropy 
the fastest 180° flip of a moment is achieved by apply- 
ing the field H perpendicular to M^. The field is turned 
on, and the moment starts to precess around H. After 
a half of the precession period it reaches the direction 
My = — Mj, at which time the field has to be switched 
off. This scheme constitutes the simplest example of 
"precessional" or "ballistic" switching. In the presence 
of anisotropy precessional switching becomes more com- 
plicated, but is still possible^ - — The key ingredients of 
the precessional switching are extremely short pulses of 
precise duration and sufficiently strong field magnitude. 
Ideally their rise and fall times should approach zero, but 
this is very hard to achieve experimentally as the pulse 
duration should be of the order of picoseconds. Preces- 
sional switching was observed using the fields created by 
pulses of synchrotron radiation^ or by a careful pulse 
shaping using femtosecond lasers ji 9 . but is not yet used in 
applications. The speed of the precessional switching in- 
creases with the field magnitude, but cannot be increased 
infinitely due to the breakdown of the ferromagnetic state 
of the sampled 

Another widely discussed class of switching scenar- 
ios is the application of the high-frequency oscillating 
field, or an RF-signal. When the frequency of the RF 
signal is close to the eigenfrequency of the magnetic 
moment, a magnetic resonance occurs and a continu- 
ous precession state is established. The amplitude of 
the precession is growing as the strength of the RF 
signal is increased. Sufficiently large signal would in 
principle be able to increase the amplitude so much 
as to switch the magnetization from "up" to "down". 
Experiments observed an RF-assisted switching 21 hap- 
pening inside the Stoner-Wohlfarth astroid. In this 
case the RF signal helps the external field to drive the 
transition. To improve the RF-assisted switching the 
schemes with variable ( "chirped" ) frequency are designed 
theoretically^ - — Their goal is to keep the frequency 
equal to the instantaneous, amplitude-dependent reso- 
nance frequency of the magnet. Various instabilities may 
prevent the purely RF switching by destroying the coher- 
ent single-domain state of a sample. However, they are 
supposed to be suppressed in the magnetic particles with 
the sizes below 10 -j- 20 m m 12 ! 25 

The phenomenon considered in the present paper is 
different from all of the above. We consider conventional 
switching by a field jump, but concentrate on the switch- 
ing time r m dependence on the jump rise time r/,. The 
function r m (r^) exhibits an unexpected minimum which 



was not discussed in the literature. As it will be shown 
in the conclusions, the minimum of r m (r^) results from 
an admixture of a precessional (ballistic) switching to the 
conventional switching, an effect which can be named a 
ballistic-assisted switching. Normally, the ballistic con- 
tribution of a small perpendicular field is quenched by 
the large anisotropy, but here it is restored by the time- 
dependence of the switching field during the rise time 
of the jump. This type of ballistic contribution can be 
observed for infinitcsimally small perpendicular fields, 
which distinguishes it is from the purely ballistic case 
where a finite field comparable to the anisotropy fields is 
required. Our analysis also shows that the phenomena of 
ballistic-assisted and RF-assisted switching are compli- 
mentary. Both can be treated on equal footing by con- 
sidering the averaging of the perturbation field torque 
during the precession cycle of the moment in the strong 
anisotropy field. 



II. MODEL 

Our treatment is based on the Landau-Lifshitz equa- 
tion, and does not take into account thermal fluctuations. 
We consider a single domain magnetic bit described by 
a magnetic moment M = Moil, where n is a unit vector. 
The magnet has an easy anisotropy axis directed along z, 
and its anisotropy energy is given by E a = —{\/2)Kn 2 z . 
This anisotropy creates two equilibrium directions of the 
moment along +z and —z. 

The switching field H = H(t)z is directed along z as 
well. It favors the +z direction for H > and — z for 
H < 0. For large enough magnitudes, \H\ > K/Mq, 
there is no equilibrium in the direction opposite to the 
field. 

When the external field is pointing along the easy axis 
direction, magnetic switching relies on the fluctuations 
near the equilibrium position. Without them a moment 
pointing exactly along +z will not be switched by any 
amount of negative applied field. Instead it will remain 
at the point of unstable equilibrium, until an initial fluc- 
tuation occurs and then grows with time. The switching 
time in this case strongly depends on the fluctuation mag- 
nitude 69, being infinite for 69 = 0i 28 ' 29 In order to model 
the required initial fluctuation we apply a small bias field 
= H±x perpendicular to the easy axis. This field is 
set to be constant in time. It creates a controlled devia- 
tion of the initial magnetization from the field direction 
and makes the problem well defined. 

The switching dynamics is described by the Landau- 
Lifshitz-Gilbert (LLG) equation. 



M = -7 



dE 
dM 



x M 



o 



-[MxM] 



where 7 is the gyromagnetic ratio, E — E a — (H+H^)-M 
is the magnetic energy, and a is the Gilbert damping 
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FIG. 1: Top: Stoner-Wohlfarth astroid in the ff-space and the 
trajectory of the magnetic field change (dashed line) during 
the rise time of the pulse with axial field sweep from +Ho to 
—Ho at a constant perpendicular bias H±. The lower panels 
(A) and (B) correspond to the initial and final directions of 
the field. They show the "up" and "down" minimum points 
M u , Md, and maximum points X of the magnetic energy. 



constant. In terms of the unit vector n it reads 

de 



<9n 



x n 



an x n 



(1) 



with e = -fE/Mo. Using the spherical an- 
gles (0,<f>), denned so that n = {n x ,n y ,n z } = 
{sin 9 cos 4>, sin 9 sin <fi, cos 9}, one gets a system of equa- 
tions 



(1 + a 2 ) 
(1 + a 2 ) 



1 de 
sin 9 d(f> 
1 de 



a 



de 
89 ' 
a de 



sin 9 09 sin 2 



(2) 
(3) 



Actual materials are characterized by a <C 1 and we will 
always calculate up to the linear terms in a, e.g., 1+a 2 ~ 
1. In our case the energy has the form 



U)Q 2 

— nt 

2 z 

L0 Q 



h{t)n z 



h i n 7 



(4) 



= — ^- cos 2 9 — h{t) cos 9 — h± sin 9 cos < 



where ojq = jK/AIq, h — jH, h± = jH±. The smallness 
of the bias field is ensured by the condition h± < w . 

The field sweep is assumed to be linear in time and 
given by the expressions 



h(t) 
h(t) 
h(t) 



+h Q , (t < 0) 



2f 



h(t) = h I 1 - — j , (0<t< T h ) 
-h , (t > Th) ■ 



(5) 



The trajectory of the magnetic field change and the posi- 
tions of the minima and maxima of the magnetic energy 
for the initial and final field orientations are shown in 
Fig. [1] We start with the positive field h = +ho > ojo 
which guarantees that the magnet is initially pointing 
close to the +z direction. This state will be called an 
"up-equilibrium" . At the end of the reversal the mo- 
ment reaches the "down-equilibrium" , corresponding to 
h = —ho- The initial and final directions of the mo- 
ment are determined from the conditions de/dcf> = 0, 
de/89 = with h = ±ho. This gives <fi — and an 
equation for 9 reading 

u>o sin 9 cos 9 ± h sin 9 — h± cos 9 = 0. 

In the limit h± <^ ujq one finds the following approxima- 
tions for the values of the polar angles in the up- and 
down-equilibria 



uj + h 



lo + h 



(6) 



As the field is swept from positive to negative values, 
the up-equilibrium disappears and the magnetic moment 
starts to move in a spiral fashion towards the down- 
equilibrium, approaching it exponentially. To define a 
finite switching time we have to introduce a provisional 
cut-off angle 9 SW and calculate the time it takes to reach 
9 SW during the switching process. The remaining dis- 
tance from 9 SW to 9 a takes extra time, but this time in- 
terval does not depend on the field sweep time since in the 
regimes studied in our paper < r m , and the remaining 
motion happens at a constant external field. We use the 
commonly adopte d 1 ' 8 ' 10 ! 12 value of 9 SW =w/2. 



III. NUMERIC RESULTS: NON-MONOTONIC 
DEPENDENCE OF THE SWITCHING TIME 

The LLG equation can be easily solved numerically and 
the switching time dependence T m (rh) can be obtained. 
Fig. [2] shows the results of such modeling for a particular 
parameter set (see figure caption). The minimum of r m 
is clearly observed. 

We will show in the subsequent sections that the 
switching time can be approximated by an expression 



3h Q + loq ln[ho/irh 2 ± T h ] 



4ho 



-Th 



2a(h Q - uj ) 



TR(a,ho,uo) 



where tr is independent of Th and h± . As one can see in 
Fig. [2] the correspondence between the actual (numeric) 
and approximate curves is quite good and reproduces the 
minimum of r m . 

A note on the numeric calculation is due here. It is 
more convenient to follow the time dependence of the to- 
tal energy e(t), than that of 9(t). The reason for that is 
as follows. In the regime considered here the switching 
time and field sweep time satisfy r m > Th and the switch- 
ing threshold 9 = ir/2 is reached when the external field 
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FIG. 2: Switching time as a function of field sweep time (time 



is measured int he units of ui n 



Parameter values are a 



0.01, ho = 3.5 ujq, h± = 0.001 ujq. The thin solid line is the 
approximate expression derived in this paper. 



is already time-independent, h = —ho- According to the 
LLG equation, at constant external field the time deriva- 
tive of energy is strictly negative 

e(t) = -ah 2 < , 

and e(t) is a strictly decreasing function of time. This 
property greatly simplifies the solution of the equation 
£( r m) = £swi where e sw is the new cut-off, introduced 
instead of 9 SW . Note that in contrast to sit) the time 
dependence of 9 it) is non-monotonic. One can easily un- 
derstand that by recalling that in the a <C 1 limit the 
spiral motion of the moment approximately follows the 
equipotential lines. Due to the presence of the bias field 
h± the latter differ from the 9 = const lines, leading to 
the oscillations of 6{t) in time. In other words, the fact 
that e(9, <f>) depends on both spherical angles makes the 
cut-offs 9 SW and e sw not completely equivalent. Never- 
theless, they serve the same purpose with the latter be- 
ing a more convenient choice. We adopt the cut-off value 
s S w = £( 7r /2,7r/2) = 0, closest to the original definition. 1 



IV. ANALYTIC APPROXIMATION FOR THE 
SWITCHING TIME 

The switching process consists of two stages. The first 
stage is the field sweep time interval, < t < Th- The 
second stage is the motion in the constant field for the 
time interval 77, < t < r m . Below we use two different 
approximations to find the magnetic dynamics in each 
stage. 



A. First stage 

The idea for the first stage approximation is to assume 
that the deviation of n from +z is small. The rational for 
that is provided by the following argument. The initial 
position of the moment is given by 9 U ~ hx/i^o + ho) <C 



1, i.e., is very close to +z. We then assume that proximity 
to the +z holds throughout the first stage if the sweep 
time is not too long. The precise condition imposed by 
this assumption on 77, is not clear at this point but will 
be obtained after we do the calculations. 

According to the above, we linearize the LLG equa- 
tion near the +z point. In the linearized equation the 
unknowns are the two projections in x ,n y ) of the unit 
vector. Both are small for n close to +z. One gets a 
linear system 



n u 



-aiu>o + h)n x — (wo + h)n y + ahx 
- (wo + h)n x — a{uj + h)n y — hx . 



Introducing a notation ui(t) = ujq + hit) we rewrite it as 

(7) 



— aui(t) —u){i) 
ui{t) —au>it) 



ah± 
-hj 



The matrix on the right hand side can be diagonalized by 
changing variables to £ = n x + in y , r\ = n x — in y . Using 
them we get two decoupled equations 

i = (i - a)[u(t)£ - h x ] , 
f) = + a) [w(t)j7 - hx] , 

which turn out to be complex conjugates of each other. 
Consequently, we can solve either one of them. Denoting 
fi = i — a, we search for the solution of the first equation 
in the form 

£(t)=Ait)e fi ti"^ ds . 

For future notation we define a phase function tp(t) = 
f*ujis)ds. The solution is found to be 

f (t) = £(0)e w(t) - nhxeP*® f e _w(u) d« . 



Going back to (n^, n y ) we obtain after the necessary al- 
gebraic transformations a solution 



n x = e av {n x o cos ip — n y o sin ip— 

- hx\iS- aC)} cos ip-iC + aS)} sin 95] } 
n y = e~ aip {n x o simp + n y o cos ip— 

- hxiiS- aC)} sin p + (C + aS)} cos p] } 

where we have defined 



(8) 



Sit) 
Cit) 



s^W S m<j>is)ds , 
3 av(s) cos <f>{a)d8 . 



(9) 



In our case the initial conditions are given by 
hx 



n x0 



, 



(10) 
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thus both projections n x and n y are proportional to h± 
and we should be able to satisfy the assumption of small 
deviation from the origin for sufficiently small bias field. 
Below we will calculate how small should h± be to ensure 
small deviations in Stage I. 

For the linear field sweep ([5]) the phase <p(t) is a 
quadratic function 



Mo + h 1 



2.s 

n 



ds = (ujq + ho)t — hp 



In this case the integrals ([9]) can be found exactly and ex- 
pressed through the error function of complex argument 
(Appendix©. 

Here we will consider a useful approximation valid 
in a large region of parameters. The phase ip(t) has 
one maximum on the interval [0, Th] at the point t m — 
Th(&o + ho)/2ho (recall that luq < ho). The presence of a 
maximum means that the integrals for C{th) and S(Th) 
can be approximated by a steepest descent (stationary 
phase) method in the case of a large change of ip on the 
integration interval [0, Th]. This certainly requires the in- 
equality u>oTh 3> 1 to hold but, as it turns out below, 
sometimes an even stronger condition is needed. The 
steepest descent calculation is performed in Appendix lAl 
and gives an approximation 



C 



KT h 

ho 

TTTh 

h 



sin [4> m - ~J + ^ cos 



7r\ a 
— J — — sin 
4/ 2 



<p(t T , 



(ujq + h ) 
4/i 



-Th 



(11) 



Substituting this into flSJ), using ^(r/j) = ujo t h and per- 
forming the calculations we find 



n x (T h ) 



n y {T h ) 



luo + ho 



cos ujot ■ 



e ^h ± 4 / ^ (sin A<p - £ cos Aip) 
ho 2 



- au "> Th — smujpt- 



luo + ho 



(12) 



-e aAv h A 
where we have defined 



TTTh -r — a —r 

—— (cos Aip + — sin Aip) , 
ho 2 



A U \ ! \ (^0 - hp) 2 

A(p = <p(t m ) - (p(r h ) = T h 



iho 



(13) 



= Aip- 



We observe that in both formulae the first term on the 
right hand side is initially small and further decreases as 
a function of Th, while the second term increases with 
Th- Therefore the condition for small deviations can be 
formulated as the smallness of the second term 



TTT h 



hp 



Explicitly separating the product LOpTh, we can write the 
condition on the bias field 



ho 



1 



itujo V uoTh 



■ exp 



(wo 



ho) 2 



Ah uo 



(uoi~h) 



which will guarantee the validity of the small deviations 
assumption. Additional inequalities (|A4[) enabling the 
approximation are listed in Appendix [A] and have to 
be satisfied as well. We will return to their discussion in 
SecM 



B. Second Stage 

During the Stage II the external magnetic field is con- 
stant, H = —Hoz + H±x. The action of the bias field 
H± can be viewed as a perturbation of the axially sym- 
metric problem with H± = and H = —Hp. For the 
unperturbed problem the switching time is a know n 28 i 29 
as a function of 8i n , the angular deviation of the magne- 
tization from the easy axis at the beginning of Stage II. 

To find the perturbation corrections to the axially sym- 
metric problem we employ the method of deriving an ap- 
proximate differential equation for the total energy s in 
the limit of small Gilbert damping constant (see, e.g., 
Ref . |30| ) . In the a — > limit the motion of the moment 
M(t) can be viewed as a fast precession along the e = 
const lines and a slow motion from one equipotential or- 
bit to the next one nearby. Up to the linear terms in a 
the change of energy upon one precession cycle around 
an orbit is given by 



As 



de 
On 



dn 



-afi{e) 



and the period of this cycle is 
Jr \de/dn\ 



= h(e) , 



(15) 



(16) 



where the integrals are taken along the constant energy 
orbit r(e) on the unit sphere. In this approximation the 
differential equation for e(t) reads 30 



de 
~dt 



Ae fi(e) 



(17) 



Its solution is given by 



de 



and determines the time required to move from the orbit 
with energy e = e\ to the orbit with e = £2- 

To find the integrals (fT5)) and ([To 7 ]) , we have to find the 
orbits r(e). When h± is small, one can expect that the 
equipotential orbits will be close to those in the unper- 
turbed case with h± = 0. The latter are the circles of 
constant polar angle. Indeed, in the absence of the bias 
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and 



FIG. 3: Rotated reference frame and the definition of 8*. 
Unperturbed equipotential orbits are shown in dashed lines 
and the actual orbits given by the solid lines. Grey arrows 
show the transformation of the representative orbits as the 
bias field h± is turned on. Point X is the energy maximum 
position in the presence of bias field. 



field e = Sq(9) and for any given energy the polar angle 
is given by an inverse function 9 = 9q(s). This statement 
is true most of the time, however important exceptions 
exist. As shown in Fig. [3l the orbits are indeed relatively 
close near the equator. But near the North pole N of 
the sphere the orbits of the perturbed energy are small 
circles around the maximum point X, while the original 
orbits are small circles around N. They are not relatively 
close in the sense that the perturbation is larger than the 
orbit size. 

One can remedy the situation by switching to a coordi- 
nate system with the z' axis going through the maximum 
point X. The new coordinate system (a/, y' , z') is rotated 
with respect to the original (x, y, z) system by an angle 
6* around the fixed y 1 — y axis (see Fig. [3]). The spher- 
ical angles of the new system will be denoted by 9' and 
(f>' . The advantage of the rotated system comes form the 
fact that the perturbed orbits are close to 9' = const 
circles everywhere, except in the vicinity of 9' = n. 
For the region of interest < 9' < ir/2 it is possible 
to define the shape of the orbit r(e) as a perturbation 
9'{<j)' ,e) = 9o(s) + h±S9((j)',e) around the circles of con- 
stant 9' . 

To prove the geometrically intuitive statement of the 
preceding paragraph one has to invert the equation 
e{9',<f)') = e. The form of energy function in the new 
spherical angles if calculated in the Appendix [Bj Up to 
the first order in h i we find 



e = e Q (9')+Pe 1 (9',cf>') + 
where the small parameter is 



h - uj 



(18) 



(19) 



£o(#') = -y cos 2 <9' + /i cos6i' 
£ X {9') = w o (l-cos6>')sin6>'cos0' 



(20) 



As expected, the zeroth order term is given by the un- 
perturbed energy as a function of the new polar angle 
9'. The same would have happened in the original co- 
ordinates, but there is an important difference between 
the new and original coordinates, which manifests itself 
in the behavior of E\ at small values of 9' . Note that the 
second term in (fT8|) grows with 9' slower than the first 
one: E\ ~ 0' , while £o — £o(0) ~ 6' . As a result, the 
inequality (3e\ <C £o — £o(0) holds uniformly in 9'. We 
will see in a moment that it is precisely this uniformity 
that is important. In the original coordinates as 9 — » at 
an arbitrarily small but fixed f3 the second term exceeds 
the first one, and the uniformity is violated. 

The orbit equation 9' = 6' (</)', e) is found form the 
constant energy condition 

£ O (6>')+/3(l-cos(9')sm6>'cos0' = £ , 

which has to be solved for 9' . The solution is searched in 
the form of a power series in j3 

9' = 9 (e)+(39 1 (0', e) + . . . 

where 9o(s) is the inverse function of £o(#), as was already 
discussed above. Up to the first order in (3 we find 



%{e) - P 



£i(flo(g),0') 
{deo / d6)\ 9= g ( e ) 



(21) 



In the "dangerous" limit of small 9 (e) near the North 
pole the second term in (f2~Tj) is proportional to 9o 2 (e) thus 
being a small correction. Due to j3 <C 1 it remains a small 
correction for the energy values up to near the equator 
£ « 0. This algebraically proves the geometrically in- 
tuitive conclusion made above: the perturbed orbits are 
close to the 9' = const lines. 

Using the approximations (fT8|) and (12~T1) for the energy 
and orbit shape, the integrals (|T5l) and (fTTH) can be eval- 
uated up to the first order in /3. The details are given in 
the Appendix [Cl Substituting the results into Eq. (fT7|) 
we get a differential equation 



de 
~dt 



de 
89 



(22) 



e=e (e) 



Up to the first order terms in /3 this is the same equa- 
tion as one would have for magnetic switching in the 
unperturbed case with h±_ = 0. We conclude that in 
the rotated coordinates one can approximate the switch- 
ing time by the expression for the unperturbed case. 
The latte r ' is reviewed in Appendix [D] and gives the 
switching time as a function of the starting angle 9i n 
and the cut-off angle 9 SW . To find the time ti spent 
by the moment in Stage II we just have to set 9i n to 
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be the value of 9' at the end of Stage I, and 8 SW to 
be the value of 9' at the selected switching moment 
given by e — 0. Thus 0j„ = 9'(n x (Th),n y (Th)) and 
9 S w = 9' (9 = tt/2,0 = tt/2) = tt/2. The substitution 
is performed in Appendix [D] and gives 



T'2 



1 

2^ 



In 



hp - ojq cos 9 in 

h - w "' \h Q (l - cos9 m ) 

1 / h - w cos 6> m 
■ In 



ft-o+^o \^o(l + cos6» in ) 



(23) 



C. Total switching time 

The total switching r m time is given by the sum of 
the contributions from Stages I and II and equals r m = 
Th+T2- To use the expression (|23|) for T2 we need the value 
of 9i n . This angle is given by the distance between the 
endpoint of Stage I {n x (Th),n y (Th)} and the position of 
the energy maximum point X given by (— h±/(ho— wo), 0) 
on a unit sphere. Since the point {n x (rh), n y (Th)} and 
the point X are both close to +z, we can approximately 
write 

in = \j (n x (T h ) + h ± /{h - w )) 2 + n y (r h ) 2 (24) 

A long expression for 9i n can be obtained by substituting 
the formulae (fT2"]) into the equation above. 

As it was already discussed, the expressions (fT2"|) rep- 
resent n x and n y as a sum of two terms, where the first 
decreases and the second increases with time. Calcu- 
lations are substantially simplified when the increasing 
term dominates, which turns out to be true in a large 
part of the parameter space. Even for small values of a 
this is guaranteed for l/(wo + ho) <^ \J r^/ho or 



u h 



(w + h ) 2 



(25) 



Since we already assumed that u>o T h 1, and the right 
hand side of ([23)1 is always less than 1 /4, this inequality is 
automatically satisfied whenever we can use the steepest 
descent approximation (|T2|) for Stage I. 

Next, we assume that one can also ignore h±/(ho — uJo) 
in the first term of Eq. (|24[) compared with n x (rh) given 
by the dominant term of (|T2|) . This condition is satisfied 
for 1/ (h - w ) < \frhJho or 



u T h > 



uj h 



(h - w ) 2 



(26) 



As a result, leaving only the dominant terms we obtain 



''in ~ hj_\l — e v 
h 



(27) 



Using the smallness of 9i„ , we approximate sin 9i n ss 9 L . 



and cos $i n ~ 1 in the expression (|23[) and rewrite t-i as 



T-2 



1 



2(h - w ) 



2a \ h Q - w M I, h 9 2 
1 , / ho - w 



-, In 

h + w 

Mi/C) 

2a(h - w ) 



2h 



(28) 



where 

r R (a,ho,uo) 



1 

2^ 



1 



■In 



2(h - w ) 



ho ~ w \ ^0 
1 , fh Q -ujo 



ho 



In 



2/ln 



(29) 



is a time interval independent of Th and ft.^. 

Substituting 6* in from Eq. (|2"7l into Eq. (f2"5|) . we pro- 
duce the first principal result of our paper, a formula for 
the switching time 



7"m — 



h\[ho/Trh\T h ] - aA<j)(T h ) 
2a(h - w ) 



tr 



(30) 



3/i +w lnf/io/Tr/i^r^] 



4fco 



-Th 



2a(h - w ) 



T R (a, h Q ,ui ) 



The obtained r m (r^) dependence indeed has a mini- 
mum. It is reached at the optimal field sweep time 



1 



Aho 



2a(h — w ) 3h + w 



(31) 



that is independent of the bias field. This formula is 
our second main result. The independence of of h± 
is a result of the logarithmic dependence in the second 
term of (|30[) and ultimately stems from the logarithmic 
dependence of T2 on the initial deviation angle. 

The minimal switching time T m (r^) corresponding to 
the optimal field sweep time equals to 



Tm(Th) 



1 + ln 



( a(feo-^o)(3feo+^o) ^ 

2a(h - w ) 



tr ■ 



It does depend on h±, which is a quite natural since the 
initial deviation from the easy axis is controlled by the 
bias field. 

Expression (|30[) has its limits of applicability discussed 
in the next section. In particular, it is not applicable for 
small Th where the steepest descent approximation (fTTj) is 
invalid. Nevertheless, one can easily calculate r m (0) since 
in this case there is no motion in Stage I, and, according 
to Figs. [T]and|3l the initial angle for Stage II is simply 

h± h±_ 



" Lit U — t 1 i 

h + w h - w 
Using this value of 9 in in (f2"5|) we get 



,(0) 



1 



a(h 



■In 



o - w ) 



hi - w, 



2hoh 



TR 



(32) 
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The drop of switching time from 77, 
value is given by a formula 



to the minimal 



In 



,(0) 



2ah^(3h +u>a) 



2a(h - lj ) 



Note that this difference is again independent of the bias 
field h±, as long as the approximation (f30]) is valid. The 
fractional change (r m (0) — 7"m(T^))/r m (0) will depend on 
the bias field, being an increasing function of hj_. 



V. VALIDITY REGIONS OF THE ANALYTIC 
APPROXIMATION 

A number of approximations were made in our deriva- 
tion and the final expressions can only be used in the 
region on their validity. 

The approximations made in our treatment of Stage I 
are as follows. 

(a) The steepest descent method employed to evalu- 
ate the integrals © has to be sufficiently accurate. The 
required conditions (see (|A4I) in Appendix lAl read 



> 



1 



ho ho — ojq 



> 



1 



(33) 



where the second inequality holds automatically, 
(b) Next, we assumed that the inequality 



h± 



h + uj 



< h A 



ho 



is satisfied, allowing one to neglect the first terms on the 
right hand sides of Eqs. (|12[) . However, this requirement 
is not new because it is already contained in (|33l) . 

(c) In order to make approximations in Eq. (|24[) we 
required the inequality (|25[) to hold. This inequality is 
also contained in (|3"3"j) and does not add new conditions. 

(d) Finally, the inequality ([HI should be satisfied to 
ensure small deviation of n from +z. 

Overall, the inequalities 



ho 



< hi 



< e 



-aAip(rh) 



< 1 



u>o 



(34) 



summarize the requirements for Stage I. For our treat- 
ment of Stage II we assumed the following. 

(e) The Gilbert damping constant should be small, 
a <C 1, to be able to use the orbit averaged equation 
of motion (IT71) . 

(f) The parameter describing the orbit deformation 
pTjl should be small, (3 — h±/(ho — LUo) "C 1. But this in- 
equality follows from Q34p and thus brings no additional 
restrictions. 

The requirements (f3"4"|) discussed above can be equiv- 
alently presented as conditions on 77, that have to be 
satisfied at fixed bias field h±. In this form they read 



ho 



(h - ujo)'' 



« T h « r, ( + } 



(35) 



where is a solution of 



ITh_ e -aA<p(r h ) 

h 



1 

hi 



We can now check when does the optimal field sweep 
time lie in the region of validity of our approximation. 
Using our result (|3T|) we can write 

1 



h 2a(h - uio) 
and substitute it into the requirement (|33|) . We get 



a < 



h p - U}Q 

2h 



The right hand side of this inequality is always smaller 
than unity, thus it automatically implies a <C 1. 

We also have to satisfy the condition (|14p at t)> = . 
For the exponent aA(f> one can write 



(ho - w ) 
a Ti 



2 H~^<1«1, 



ih n 8h 
and thus condition (TT4")) simplifies to 



The above inequalities on a and h± can be combined into 
a single requirement 



hi 



< a < 



h - ojo 



(36) 



h (ho - ujo) 2h 

If the inequalities p6l) are satisfied, the optimal sweep 
time t£ occurs inside the interval (1351) and can be calcu- 
lated using the formula (l3~Tj) . 



VI. COMPARISON OF ANALYTIC AND 
NUMERIC RESULTS 

Comparisons of the analytic approximation and exact 
numeric results were performed in Ref. and shown a 
very good agreement between the two. 

First, we compared the numerically calculated switch- 
ing times with the expression (|30[) . When the inequalities 
Q36p were well satisfied, the quality of approximation was 
very good (Fig. @|. As one approached the limits of the 
approximation's validity by, e.g., increasing a, the errors 
grew larger. 

Second, the numeric results for optimal field sweep 
time were compared with the formula pip . The cor- 
respondence was generally good (Fig. [SJ , although some 
visible deviations existed. They were attributed to the 
fact that the accuracy of the determination of is low- 
ered by a flat shape of the T m (rh) curve minimum. Be- 
cause of the shallow minimum, small errors in r m produce 
much larger errors in r^. 
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FIG. 4: Dependencies T m (jh) calculated for ho = 2.2 ujq, 
h± = 0.001 Wo, and variable a indicated on each panel. As a 
increases, the theoretical fit gets poorer due to the violation 
of the strong inequality PrT)) . 
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lines) dependencies of the optimal field sweep time on the 
system parameters. (A) fixed a and h± , (B) fixed ho /wo and 
h±, (C) fixed ho/ojo and a. When not varied, the parameter 
values are ho = 2.2wo, a = 0.01, h± = 0.005. 



sions (|A3[) for the integrals C and S and the formula 
for 9i n , one can obtain an excellent agreement between 
the two-stage theory and the no-approximation numeric 
simulations in this limit. 



In general, the analytic expression approximated the 
T m (Th) dependence up to a 10% accuracy in a surprisingly 
wide range of parameters. Such accuracy is certainly 
sufficient for the estimates related to the device design. 

When Th is outside of the validity region of the results 
(l3Cfl) and (|3Tj) , exact expressions ([9]) and (|24|) for the inte- 
grals and the initial angle 9i n can be used. As long as the 
deviation from the +z direction in Stage I remains small, 
they provide a good approximation for r m . Consider for 
example the case of small tv Approximation (f30f does 
not work for t>, — > predicting an infinite increase of r m , 
while the actual limit r m (0) is finite and given by for- 
mula (J32). It was checked that using the exact expres- 



VII. PHYSICAL PICTURE OF THE 
BALLISTIC-ASSISTED SWITCHING 

We now discuss the physical reason for the minimum of 
the function r m (r^) which was identified in Ref. 2 as the 
contribution of ballistic switching. We start by noticing 
that the bias field has two roles in the switching process. 
First, it provides the initial deviation from the easy axis. 
Second, it alters the equations of motion for n(t). The 
first role of h± manifests itself in our formulae in two 
ways: by providing the first terms in expressions (|12[) and 
by introducing the term h±/(ho — uq) into the formula 
(|24p . In our derivation of the switching time expression 
([50)1 we have found that both contributions are negligible. 
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FIG. 6: Average ballistic contribution of the bias field. Vector 
n(t) orbits around a parallel circle C on a unit sphere. The 
torque due to h± pushes n along the meridians of the sphere. 
In constant switching field the torque contributions from the 
diametrically opposed elements dh and dl-z cancel each other. 
For variable H(t) the cancelation does not happen (see text). 
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FIG. 7: Numerically calculated witching time with pulsed bias 
field. The parameters are set to ho = 2.2ojo, h± = O.OOIojo, 
a = 0.01 (cf. Fig. 3]). Field sweep time is fixed at = 50/wo, 
close to the optimal sweep time = 48.2/o;o. The time 
dependence of the pulsed bias field h±(t) is shown in the inset. 
As the width of the pulse At approaches the theoretical target 
of 2^/rh/ho ~ 10/wo , the switching time approaches the value 
obtained at h± = const (shown by a horizontal line). 



Therefore within our approximation only the second role 
of the bias field is important. 

Let us proceed by discussing this role qualitatively. Re- 
call that in the absence of anisotropy and other fields the 
torque x A/on due to the bias field would rotate the 
unit vector n from +z to — z along a meridian of the unit 
sphere (dashed line in Fig. [()]), in a "ballistic" or "preces- 
sional" fashion. In our case a weak bias field is applied 
on top of the strong uniaxial anisotropy and switching 
field, which induce a fast orbital motion of vector n(i) 
along the parallel circles (line C in Fig. [5J). The bias field 
still attempts to move n along the meridians, but now 
its action has to be averaged over the period of orbital 
motion. As illustrated in Fig. [(J in constant fields H||z 
averaging gives zero due to the cancelation of the con- 
tributions from the diametrically opposed infinitesimal 
intervals dl\ and dli of equal lengths. This way ballistic 
contribution of the bias field is quenched. However, the 
contribution of Hi docs not average to zero for a vari- 
able switching field H(t). In this case the velocity of n 
changes along the orbit, the times spent in the intervals 
dl\ and dZ 2 ar e different, and the contributions of the two 
do not cancel each other. We conclude that in the pres- 
ence of a time dependent field H(i) ballistic contribution 
of the perpendicular bias field is partially recovered. 

Based on the discussion above, we may expect to find 
that the largest contribution of the ballistic switching 
will happen when the change of velocity is biggest. The 
measure of the velocity change is Auj/lj, where Aw is the 
change of the velocity as one precesses from the interval 
dl\ to the interval dli- We may estimate 

Aw wT Co 

UJ LU US 2 02 ' 

where T — 2-k/lo is the instantaneous period. We see 



that the velocity change estimate diverges at the point 
= 0. But this is exactly the stationary phase point t m 
that gives the largest contribution to the integrals (O in 
our approximation. We have thus established a one-to- 
one correspondence between the qualitative description 
of the ballistic contribution and our derivation of the an- 
alytic approximation (|30|) . To illustrate the fact that 
the action of hj_ is only important near the stationary 
phase point t m , we have performed a numeric experi- 
ment with pulsed bias field h±(t) changing in time as 
shown in the inset of Fig. [7] It is kept constant for t < 0, 
switched off at t = 0, and then switched on again for a 
short interval of time At, centered around the t m point. 
With the pulsed bias field the ballistic contribution is 
only present during the interval At. It follows from the 
steepest descent calculation of Appendi x lAl tha t formulae 
(fTTj) would be valid already for At > 2y/Tf l /ho, and when 
this inequality is satisfied our theory would predict the 
same switching time (|30[) for pulsed and constant bias 
fields. The results of the numeric experiment (Fig. [7]) are 
completely consistent with this prediction. As the pulse 
width approaches the value of 2y // Tf l /ho, the switching 
time drops to the value obtained earlier at constant h±. 

The minimum of the T m (r^) function can be under- 
stood as follows. Ballistic contribution helps to move 
vector n from +z to — z and is thus responsible for the 
initial decrease of r m . As the sweep time grows larger, 
the change of the orbital velocity during the precession 
period decreases and the ballistic contribution averages 
out progressively better. The helping effect of ballistic 
switching is lost and r m starts to increase as it normally 
would. 

Finally, we want to remark that ballistic contribution 
to switching can be also viewed as a phenomenon com- 
plimentary to the magnetic resonance and RF-assisted 
switching<2i~— In the case of RF-field application the ex- 
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ternal field H is constant, while the bias held Hj_(t) is 
time-dependent. Here the contributions of the bias held 
torque on the intervals dl\ and dfo in Fig.[6]do not cancel 
each other because the torque itself changes with time. 
This, again, leads to a nonzero average of the bias held 
contribution on a precession orbit and creates a help- 
ing effect for the magnetic switching process. From this 
point of view the magnetic resonance and the time de- 
pendence of the axial switching held are two different 
ways to achieve the same goal: a non-vanishing average 
contribution of the bias held torque on an orbit. 



VIII. CONCLUSIONS 

We have identified and investigated the phenomenon 
of ballistic contribution to the conventional magnetic 
switching by a time-dependent held. An analytic approx- 
imation is derived for the ballistic-assisted switching time 
in a constant perpendicular bias field. It is also shown 
that, if practical, a constant bias can be substituted by 
short pulse of bias field applied near the stationary phase 
time point. Our results provide a convenient approxima- 
tion for the optimal held sweep time, an important pa- 
rameter in the device design. The expressions obtained 
in this study can be used as a starting point for the inves- 
tigations of the switching time in granular media, where 
each grain can be modeled by a single moment and bias 
held is produced by the other grains or by the spread of 
grain orientations. 
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Appendix A: Steepest descent approximation for 

C(t) and S(t) 



Here we evaluate the integrals ((9|) 



S(Th) 



e av(t) smip(t)dt 



C(t h ) = / e^W cos ip{t)dt 
Jo 

with the phase ip(t) given by a real quadratic function 



Lp = (uj + hp)t - hp — . 

Th 

The integrals in question can be obtained from the real 
and imaginary parts of a complex integral 



I = r e-wWdt = C-iS 
Jo 

(j, = i — a . 



(Al) 



By completing the square we can rewrite 

Qo + q;o) 2 h Q 2 

<PW = 77 T h V t m ) , 

where t m = ThioJp + ho)/2ho is the point of maximum 
phase, < t m < Th- Now 



' J 7 



I = e 4ft o 
Jo 



(A2) 



Changing variables to z = y JjJio/Th(t—t m ) , we can write 
down J as 



J 



Th 



dz 



r 



where T is a straight line in the complex plane go- 
ing between the points z\ = —\Jnho/Th t m and z-i = 
a/ [ih^jTh {jh — tm). Due to our dehnition of z, line T 
crosses the complex zero point. 

Since the integrand of J is a regular function, inte- 
gration can be performed along any contour connect- 
ing z\ and Z2- The integral can be expressed in terms 
of the error function of complex variable Erf (2) 
(2/y / 7r) exp(— z 2 )dz a: 

J-- 



,31 



-^-^(Erf^+Eri^)) (A3) 

The above in an exact formula. The steepest descent 
approximation corresponds to the case of large absolute 
values \zi \ ^> 1, \z2\ ^> 1- Due to the smallness of a one 
has \n\ ~ 1, so these conditions translate to 



—t m > 1 , 
Th 



— (Th - U, 
Th 



» 1 



or equivalently 

(hp + u) ) 2 
Ah 



T h > 1 



(h - up) 1 
4/i 



T h > 1 



(A4) 



Using 51 Erf(z) — > 1 for \z\ — > 00 we find the approxima- 
tion 



J 



fiho V hp V 2 



1 



where we have also expanded in small a. Substituting 
this back into (IA2[) we get 



KTh p j( 7r /4-y m ) 

hp 



lot 

T 



(A5) 



where ip m = cp(t m ). The real and imaginary parts of I 
give C and S according to Eq. cnt 



C 

s 



hp 
KT h 

hp 



i(«^m-|)-|sin (4>m~ J) 
($m ~ ^) + I COS (j) m ~ ~) 
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Appendix B: Energy in rotated coordinates 

The relationship between the projections of n and h in 
the primed and original coordinate system are given as 



-n z sin (/* + n x cos tf* = 

- cos 9' sin 6L + sin 9' cos 6' cos ( 



and 



h' z = h cos 9* — h± sin 
h' x = h sin + h ± cos 0* 

One can now rewrite Eq. Q through the angles (#',</>') 



£ = ~y n * -n x h x -n z h z = 
wo , 



= — ^ (cos 9' cos + sin 9' cos sin 0* ) 2 

— /i (cos 9' cos 0* + sin 6*' cos 4>' sin ) (B 1 ) 
— h± (sin 0' cos 0' cos 0* — cos 9' sin 0* ) 

The above is the exact expression. We are looking at the 
case h = —ho and h±_ — > 0. Up to the first order in h± 



sin6L 



h 



;, COS0, « 1 



Expanding the energy up to the first order in h± we get 
s^so{9')+/3e 1 {9',^) (B2) 
where the small parameter is 



h - u! 



(B3) 



and 



e {9') 



UJ 2 



cos 2 9' + h cos 9' 



st(9') = LJ (l-cos6')sm6'cQS<f>' (B4) 

The first term in the expansion (|18[) is the energy unper- 
turbed by the bias field, evaluated at the new polar angle 



Appendix C: Integrals along the perturbed orbits 

First, we calculate the approximate value of |cte/<9n| 
Taking the identity 



de 



dn 



de 



l 



sin 2 9' 



de 



and expanding in small j3 up to the first order we find 
de 



dn 



de Q dei 
39' P 89' 



(CI) 



Next, we need the element of orbit length \dn\. Using 

\dn\ = \J sin 2 9d(f)' 2 + d9' 2 and calculating up to the first 
order in /? we get 

\dn\ = (sme o -0™ 6 °^ 6o > <l, ' ) )d<P> 



(C2) 

V {deo/dd)\e=8 ( E ))^ 

Using the form of e\ (|B4|) we will rewrite it as 

\dn\ = (sin 6»o - /3A(e) cos<p') d<f>' 

We will further use a notation 

e 1 {9o{e)^') = B(e)cos4>' 

To perform the integrals (IT51) and (|16l) we use the ex- 
pansions (IC1|) and (|C2p . and, expanding up to the first 
order in /?, get 



dE 


dn = 






dn 




Jo 



de Q 
d9 



sin 9q- 



' dB deo 

P ( -Qj S[n9 " - ~qJ A ) C0S< ^' 



The integral of first order term in j3 vanishes and we get 
dE 



dn 



dn = 2vr% sin 9 + 0(/3 2 ) 
d9 



Perform a similar calculation for the integral (|16[) we get 



dn 
r \de/dn\ 



2tt 



sin 6>o 



G(/3 2 ) 



de / 86 

According to (fT7f the differential equation on e(t) reads 

2 



de 



de 
d9 



0{fi 2 ) 



(C3) 



8=8 (e) 



Appendix D: Switching time in the unperturbed 

case 

The problem of switching time of a uniaxial particle 
in the absence of perpendicular field was probably first 
solved by Kikuchi 2 ^ for the case of H z = 0. The deriva- 
tion was generalized to arbitrary H z by many authors, in 
particular in the appendix of Ref. [29|. Here we re-derive 
this result for the completeness of the presentation. In 
the case of h±_ = one can find the switching time either 
by solving Eq. (f2"2"| (same as Eq. (|C3p ) truncated to zc- 
roth order, or by a direct inspection of the system ©, 
([3]). Since £o depends only on 9, it is enough to consider 
the first equation which reads 



de 



-a(wo cos 9 + h) sin 9 
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Integrating we get 
—at — 



dd 



g sin 6{bjQ cos 9 + h) 
A variable change x = cos 9 gives 

t _ i dx 

a J Xin (1 — x 2 )(uj x + h) 



1 



In 



/ 1 



x 



1 

2a { h + wo \u>ox + h 

1 / 1 + s 
■ In 



h — ujq \0Jqx + h 



In application to our problem 9 S 
h = —ho which gives 



7r/2 (x sw = 0) a n d 



(Dl) 



1 

2a 



h 
1 



In 



hp - a; cos 9 m 

~n \/l (l - COS^n) 

.' /l - w cos in 

In 



Iiq+ujq \h (l + cos9 in ) 



(D2) 
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